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Abstract I explain the concept that light quarks diffuse in the QCD vacuum following 
the spontaneous breakdown of chiral symmetry. I exploit the striking analogy 
with disordered electrons in metals, identifying, among others, the universal 
regime described by random matrix theory, diffusive regime described by chiral 
perturbation theory and the crossover between these two domains. 
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Introduction 

In these lectures I review spectral aspects of the mechanism of the sponta¬ 
neous breakdown of the chiral symmetry in Quantum Chromodynamics. Most 
probably, the spontaneous breakdown of the chiral symmetry is a collective 
phenomenon caused by the microscopic disorder, in striking resemblance to 
the diffusive phenomena appearing in disordered metals. Despite the micro¬ 
scopic, detailed nature of this disorder is still unknown, the constraints arising 
from realizations of chiral symmetry in QCD are so strong, that allow us to 
predict several non-trivial consequences of this phenomenon. Among most 
profound, are the predictions on the spectral properties of the Dirac operator. 
As such, these are amenable to test using the lattice calculations. 

The lectures are intended to be elementary and self-contained. The outline 
of the lectures is as follows: 

• In Part 1,1 introduce the basic facts on symmetries and anomalies of the QCD, 
leading to fundamental low energy constraints. 

• In Part 2,1 explore the analogy between the metal viewed as a complex quan¬ 
tum system sharing universal properties with so-called chaotic systems from 
one side, and Euclidean QCD with light quarks diffusing in disordered medium 
built from the lumps of the gauge field, on the other side. In particular, I identify 
the relevant scales and determine the "diffusion constant" of the QCD vacuum. 


* Lectures given at the Cargese Summer School, August 6-18, 2001. 
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Then I comment on the close relation of this picture to several known descrip¬ 
tions of chiral symmetry breaking. 

• In Part 3,1 demonstrate how the hierarchy of spectral scales explains the ori¬ 
gin of the appearance of random matrix ensembles in QCD. Then I show how 
random matrix regime regime breaks down at a certain spectral scale (analog 
of Thouless energy), leading to various versions of chiral perturbation the¬ 
ory. I also point at some explicit confirmations of this picture coming from 
the recent lattice calculations. Next, I suggest how this picture can be modi¬ 
fied/generalized in the presence of external parameters of the QCD. I mention 
the possibility of (multi-critical) scaling. I conclude with an (incomplete) list 
of works on disorder and QCD, that could be used as a guide for further reading 
and also as a good starting point for the original research. 


Part 1 

Let me refresh here some low energy theorems of the QCD, which will be 
needed for the second paid of these lectures. Then I switch to the Euclidean 
regime, and I recall the Banks-Casher relation. 

QCD 

Quantum Chromodynamics is a Yang Mills theory based on the local gauge 
group SU( 3) co i or . This means that all six quarks of different flavors interact 
with gluons, as well gluons interact with themselves, with the same universal 
coupling g. Looking at the typical masses of the quarks, we see that they cover 
several orders of magnitude - u,d,s,c,b,t, ~ 5,7,150,1400,4400, 175000 MeV, 
- respectively. The origin of this hierarchy goes beyond the strong interaction 
sector and is unknown. The fact, that three masses arc less or at most comparable 
to the scale of the strong interaction, and three others arc well beyond the scale, 
suggests a simplification, in which we put three heavy masses equal to infinity, 
and three light equal to zero. The heavy sector decouples from the light, and the 
light one reveals a series of essential symmetries of the theory. We denote this 
idealization of the QCD as QCD X . First, let us note that such classical theory 
lacks any dimension-full parameter, therefore it is scale invariant. Second, if 
we introduce the notion of left and right-handed quarks, 


< 27 ? = + 75)7 QL = -(! - 75)9 


( 1 ) 


we see that in the presence of massless quarks the left and the right handed 
quarks interact independently with the gluons, i.e. we have two decoupled 
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chiral copies of the initial theory: 

Sqcd = [ dtd 3 x Y JjPWA^jP + (R^L) + S glue (A) (2) 
J f 

Constructing the standard set of the Noether currents corresponding to the sym¬ 
metries yields: 


K 

= qi^t a q 

K 

= qi^q 

K 

= qi^istaq 

K 

= qi^i^q 

tM 

scale 

= Xu e^ 


where a = 1,2,3, ...8 for SU(3)Q avor , t a arc corresponding generators (here 
Gell-Mann matrices, modulo normalization ), q = (u. d , s), C-V" 7 denotes sym¬ 
metric traceless energy momentum tensor, and we used the linear transformation 
for the original left and right currents ( V = R + L, A = R — L). 

For three light flavors, we arc left with a nonet of the conserved vector 
currents, a nonet of the conserved axial currents, and conserved currents corre¬ 
sponding to the scale invariance. 

SU{3) a x SU{3) v x U{1) A x U{l) v x [SCALE INV.] (4) 


This classical picture is however strongly distorted at the quantum level. 
First, QCD V has anomalies - i.e. certain classical symmetries arc violated at 
the quantum level. The scale-invariance is broken ( scale anomaly), 


O tM — CjM 
"'scale — 


2/?(g) 

9 


1 

4 




renorml. 


(5) 


and the appearance of beta function introduces the scale of the strong interac¬ 
tions, Aqcd ~ 200 MeV. 

Second, the singlet component of the axial current is also anomalous 

( 6 ) 

A si mi lar anomaly appeal's in QED, with Nf /2 —> 1, g — > e, and non-Abelian 
G 1 "' 1 replaced by Abelian electromagnetic F 1 "' tensor. These are the all anoma¬ 
lies of the QCD, if no external (e.g. electromagnetic) currents are added. The 
appearance of the Gell-Mann function (3 in the scale anomaly reminds about an¬ 
other feature of QCD: the interactions between quarks and gluons get stronger 
at smaller energy scales, invalidating the perturbative calculation in low en¬ 
ergy domain, contrary to the precisely opposite behavior in QED. Quarks and 
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gluons interact strongly forming colorless states, and the unraveled nature of 
long-wavelength limit of these interactions is usually coded under the name 
"confinement". 

The closer look at the experimental spectrum of elementary excitations 
(hadrons) of the QCD vacuum shows that the picture is more complicated, 
and another, on top of confinement, nonperturbative phenomenon has to take 
place in QCD. If we look at the remaining symmetries after quantizing the 
QCD X , we see that they arc 

SU{3) a x SU{3) v x U(l) v (7) 

where the last one corresponds to the baryon number conservation ( 3V 11 = IP'). 
The interaction preserves baryon number and is invariant under SU(3)y and 
SU(3)a symmetries. Since the masses of up and down quarks are 20 times 
smaller than the scale of the QCD, the lightest particles in real QCD should 
show the traces of the exact symmetry (7). This is at odds with the experiment, 
which shows that the light vector-like particles differ from axial-like particles, 
despite similar flavor content - e.g. p(770) is much lighter then axial ai(1200). 
This asymmetry holds also for baryons of opposite parity, (nucleon n(940) 
versus N( 1535)Si i) and manifests dramatically at the level of the scalars - the 
lowest pseudoscalar (pion) 7 t( 140) seems to not have a narrow chiral partner 
at all - when comparing pion to /o(400 — 1200) with full width 600 — 1000 
MeV, an interpretation of /o as a particle is controversial... This suggests 
that the vacuum state of the QCD is not respecting all the symmetries of the 
interaction of the QCD. This phenomenon is called a spontaneous breakdown 
of the symmetry. Then, chirally invariant interactions of the QCD X , acting on 
chirally non-invariant vacuum, can indeed produce such an asymmetry in the 
hadronic spectrum. 

The important hint comes from the observation by Vafa and Witten [1], that 
theories with vector-like couplings (e.g. QCD) cannot break spontaneously 
vector symmetries. We arc therefore left with the alternative: 

Q a A |0 >= | PS a >/ 0 (8) 

where Q A arc axial charges corresponding to the currents A%. Vacuum state 
respects therefore only vector symmetries. Since the Hamiltonian of the QCD X 
still commutes with all the 16 generators Qy, Q A we see that the vacuum 
state ( H |0 >= £b|0 >) is degenerated with the octet of the states Q A \0 >> 
(a = 1, 2, ...8). Indeed 

H\PS a >= HQ a A \0>=Q a A H\0> = E o Q a A \0 >= E 0 \PS a > 

(9) 

This is a basic message of the Goldstone theorem [2]. To each broken generator 
of the axial current corresponds a massless, spinless excitation corresponding 
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to quantum number of the generator. In the case of three flavors, the Goldstone 
theorem predicts the appearance of an octet of massless, pseudoscalar mesons. 
They correspond to massless pions (isotriplet), kaons (two isodoublets) and an 
eta (isosinglet). 

A priori, from the point of view of chiral symmetry, QCD allows (at least) 
two phases: asymmetric (Nambu-Goldstone phase), described above, and the 
symmetric one (Wigner-Weyl phase), where vacuum is respecting all symme¬ 
tries, Qy |0 >= Q a A |0 >= 0. Hence, a phase transition may happen between 
these two phases. Such phase transition can be characterized by the appearance 
of an order parameter. The lowest-dimensional order parameter is the expecta¬ 
tion value of qnqr + (join- It carries zero baryon number, is a scalar (vacuum 
respects space reflection), and is diagonal in flavor. Note that e.g. < ud> is not 
invariant under rotations generated by isospin matrices, therefore is not invari¬ 
ant under the vector, unbroken group of the vacuum. There are infinitely many 
other operators, which arc order parameters, e.g. (IlQrQrQl, but they carry 
higher canonical dimension. The lowest dimensional order parameter <qq> 
is called "quark condensate". It is believed that under the action of external 
parameters, like e.g. temperature and/or density, the crossover to other phases 
of QCD is possible (see lectures by Rob Pisarski and Krishna Rajagopal), in¬ 
cluding the phase with restored (approximate) chiral symmetry. It is important 
to mention, that despite some experimental signals point at the restoration of the 
chiral symmetry (cf.[3]). there is, in my opinion, no "smoking gun" evidence 
for this phenomenon and most of our understanding of chiral phase transition 
comes from the lattice studies (see lectures by Frithjof Karsch). We should also 
remember, that in real QCD, the masses of u,d,s quarks are light, but non zero, 
therefore on top of the phenomenon of spontaneous breakdown of the sym¬ 
metry we have also an explicit, albeit small, explicit breakdown of the chiral 
symmetry due to the explicit presence of mass terms mqq in the Hamiltonian. 
Therefore the chiral restoration is not really a phase transition, but rather a 
crossover process. 

Condensate, pion, GMOR 

We still consider a QCD r Let us come back to the Goldstone theorem, 
and show that indeed the presence of the condensate forces the presence of 
massless excitations in the spectrum. To demonstrate this, we use the Ward 
identity, relating the correlator of the axial SU(3) current ,4(' and pseudoscalar 
SU(3) current P b = iqt b '^q to scalar densities, 

d ** < 0|TA^(.x)P j (0)|0 >= -id(x) < 0|g[f,^]+q|0 > (10) 

Exercise . Justify (10), making use of equal-time commutation relations ( cf the 
lecture of Jean-Paul Blaizot). 
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Since Lorentz invariance implies that Fourier transformation of the correlator 
has to take a form 

J dxe ipx < 0|T^(x)P J '(0)|0 >= p^'Hp 2 ) (11) 

the Ward identity says that jr 1 V 1 (p 2 ) = —d ,J < qq >. (Taking a derivative 
coiTesponds to multiplying the matrix element by ip Note that the vacuum 
expectation value of the second term in VP = 6^ + d}i k t k has to vanish due to 
the invariance under vector transformations). Since the condensate (expectation 
value of r.h.s. of (10)) is non-zero in the broken phase, the spectral function 
n(p 2 ) has to contain a massless particle corresponding to the pole at p 2 = 0. 

Now we can make use of a powerful condition, (formulated by’t Hooft), 
relating fundamental theories with theories in which particles arc bound states 
of the fundamental constituents. This condition, known as anomaly matching 
condition, states that a composite particle has to reproduce exactly the anomaly 
present in the fundamental theory. 

To see how this condition works let us consider the axial SU(S)a cur¬ 
rent corresponding to 7r°, i.e. A 3 = g7/A 3 75(/, with q = (u, d, s ) and f 3 = 
diag(l, —1,0). This current is anomaly free in QCD, but if we allow quarks 
couple to photons, the matrix element 

< 01^^177 > = < 0|<9 At -u7 At 75'u|77 > - < 0|75(^77 > 

= |^ c (4 - e 2 d ) < 0| FF | 77 > (12) 

It is easy to understand the r.h.s. of (12). This is simply the difference of two 
electromagnetic anomalous U(l) currents, corresponding to the charges with 
fractions e u = 2/3 and e f j = —1/3, respectively. This equation represents the 
contribution from the celebrated triangle anomaly. 

From the above identity we infer immediately the form of the matrix element 

T)i^ cy 00 ~ 

< 0| A 3 177 >= -i-^—N c [e\ - e 2 d ) < 0 |FF | 7 7 > +... (13) 

where p is the momentum of both photons and the dots stand for less singular 
terms. The pole here is simply the remnant of the massless quark circulating in 
a triangle graph. If we now assume confinement, anomaly matching condition 
says that in the hadronic (composite) world, there must exist a massless, color¬ 
less object which couples to the axial current and photons and which reproduces 
exactly the pole in the anomaly above. Since massless baryons arc excluded 
(cf. Vafa-Witten theorem), this object has to be a meson. Mesons fulfill the 
anomaly condition in an easy way. We do not need to circulate them inside the 
triangle loops like quarks, it is enough if there is a relation between the axial 
current and the pion. This well-known relation reads 

< 0|A“|7T b >= i5 ab F nPfi 


(14) 
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where F\ is pion decay constant. So anomaly matching means, that we read out 
the r.h.s. of (13) as a product of (14), massless pion propagator — i/p 2 and the 
amplitude for ir ° —> 77 decay (rest of the r.h.s. of (13). The similar matching 
of the trace anomaly between the quark and hadronic worlds was discussed 
in Dima Kharzeev lectures). Let us come back to Ward identity (10). In the 
p 2 —> 0 limit the correlator is fully determined by the condensate, 

J dxe ipx < 0|TAj,(x)pJ(0)|0 >= ~^<qq> S ij (15) 

Saturating the correlator with pion states, and using (14), we read out comparing 
left and right h.s, the value of the another matrix element 

< 0\P j \7T k >= 5 jk < f > (16) 

-^7T 

With the help of the above relations we can now prove the Gell-Mann Oakes 
Renner relation [4]. Introducing the small quark masses, we see that Dirac 
equation implies 

d^{u{x)^^d{x)) = (m u + m d )u(x)i'y 5 d(x) (17) 

Calculating the vacuum to pion matrix element of the above identity , and using 
formulae (14) and (16), we get 

M 2 F 2 = (rn u + rn d ) < qq > (18) 

This relation shows how the gap in the Goldstone mode appeal's due to the 
presence of the nonzero quark mass. 

Euclidean world 

In this subsection we present the formulae that allow the transcription from 
Minkowski to Euclidean space. The advantages of working in Euclidean space 
are two-fold: first, several mathematical operations are well defined, second, 
the formulation is comparable to the lattice simulations. The following set of 
rules defines the transition (l.h.s. denotes Minkowski, r.h.s. denote Euclidean) 


■ Space time: ix 0 = X 4 , Xj = 27 


■ Vector potentials: Aq = T 44 , A, = A % 

■ Gamma matrices: 70 = 74, 7 i = 75 = 75 

■ Fermi fields: iq = q, q = q 

■ Action: iS = —S 



Then, Euclidean Dirac matrices obey 7 ^ 7 ,, + 7 ^ 7 ^ = 2() /J|y , all four matrices 
arc hermitian, as well as the Dirac operator i( 9 /t 7 /( . Finally, the Euclidean action 
for QCD reads: 


S 


d A x 


_ /lO /7(t 
^ ^ fllZ^ fllZ 


q(i'y tl D^ l (A) + im)q 


(19) 


and partition function Z = exp (-S). 


Banks Casher relation 


We will now demonstrate that the chiral condensate is related to the Dirac 
operator spectrum in Euclidean space-time. Consider a Dirac propagator in 
a Euclidan box V = L 4 in the presence of the gluonic background field ,4,,. 
From the action above, we read that propagator is Sf = (Jp(A) + m) -1 . We 
can write then: 

^<0|?(0)g(0);A|0> = -TR S F (x, x) = - ^ — ( ^ - - 

—» lim i [ dXp(X)- - (20) 

V—>00 J X — im 


where TR denotes trace over coordinates, color, spin and Dirac indices, and 
spectral density is defined as p( A) = J2k ~ ^fc) where / I/Xp. = Xkqk■ Now, 
we average the above equation over the gluonic configuration weighted with 
the full QCD measure (including standard gauge fixing etc.) 

< (...) >qcd= J[dA\(...)det(]f)(A) +m)e~ s s lue (21) 

As a result, we get 

V <qq>= i [ dX ^ X } (22) 

J X — im 


where p{X) =< p >qcd is the average spectral density, i.e. the density p(X) 
averaged over the full QCD measure. As a final step, we take a chiral limit and 
use the relation 


lim - 

»o A ± im 


= PV— in5(X) 
A 


(23) 


As a result, we relate the chiral condensate to the average spectral density around 
zero eigenvalues 


V < qq >= — 7 rp( 0 ) 


(24) 


Note that the contribution from principle value (PV) paid drops, since due to 
the chiral property [Ip, 75]+ the eigenvalues come in pairs (corresponding to 
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eigenfunctions and 75 %), so the average spectral density is an even function 
p( A) = p{— A). Note that the same property guarantees the positivity of the 
QCD measure, despite at the first look JpA-m seems to create complex measure. 

The relation ( 24) was first suggested by Ranks and Casher [5]. (For the 
careful discussion of the UV part of the spectrum, see [ 6 ]). It is very important 
that the chiral limit m —► 0 is taken after the thermodynamical limit V = 
L 4 —> 00, for otherwise the average spectral density would be zero. The result 
states that in a vector like theory with chirally symmetric spectra, the quark 
condensate is related to mean spectral density at zero virtuality (i.e. at A = 0). 

We can now ask the crucial question, what kind of mechanism can cause 
spectral density to be non-vanishing at zero virtuality? Note that the levels of 
free particle closed in the box scale like 2irn/L, with n integer, so the mean 
level spacing goes like I A 1 , and the average spectral density (proportional to 
the inverse of the mean level spacing) for free particle scales like L = V 1 ^ 4 , 
therefore will never be able to balance in the thermodynamical limit the l.h.s. of 
(24) provided the condensate is non zero. The only solution is that spontaneous 
symmetry breakdown requires enormous accumulation of the eigenvalues in the 
vicinity of zero, with the condensation rate scaling like V, so p( 0) ~ V. This 
at the first look obvious fact was first emphasized and exploited by Leutwyler 
and Smilga [ 6 ] and forms the cornerstone of the spectral analysis of the Dirac 
operator. 

We will devote the next lecture to unravel the most plausible microscopic 
mechanism responsible for such a spectral behavior. 

Part 2 

In this part, we outline the basic concepts of diffusion and translate them 
into the Euclidean QCD language, identifying in this way the hierarchy of the 
spectral scales of the Dirac operator. 

Primer on the diffusion 

The Banks-Casher relation is reminiscent of the Einstein relation describing 
the conductivity a of degenerate gas of electrons, 

<7 = 2e 2 Dp(Ep) (25) 

where D is a diffusion constant, p(Ep) is the density of states (per spin direc¬ 
tion) at the Fermi surface. Is it possible that the spontaneous breakdown of the 
chiral symmetry has also the diffusive nature, with the Fermi surface replaced 
by the zero virtuality band? We will show explicitly in this chapter, that this 
indeed is the case. 

In order to prove this conjecture, we have to remind some basic facts on the 
diffusion [7, 8]. The diffusion is a process, in which a typical distance covered 
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by the diffusing particle in a time t varies as 

r 2 (t ) = Dt (26) 

where D is the diffusion constant characterizing the medium. If we consider the 
diffusive motion in a cube L d of the linear size L, it is natural to define the time 
scale t c = L 2 /D, characterizing the time during which the particle can probe the 
whole system. The energy scale corresponding to this time, known as a Thouless 
energy [9] (in units where h = 1) is E c = \/t c = D/L 2 . For times shorter than 
t c , the diffusing particle can probe only paid of the volume. For even shorter 
times (shorter then the time of mean free path between the "dirt" causing the 
diffusion, t e ) the diffusion concept is meaningless. On the other side, for very 
large times (tp ~ 1/A), the diffusive particle will always leave the volume (see 
Edmond Iancu’s lectures). When such times staid corresponding to the inverse of 
mean quantum spacing of the quantum mechanical levels, the classical concept 
of the diffusion also becomes meaningless. The above hierarchy of scales could 
be pictured by a cartoon, where we introduced also the names of four different 
regimes using the terminology borrowed from the mesoscopic physics [7, 8]: 
Quantum( t > tp), Ergodic (universal) (tn<t< t c ). Diffusive {t c < t < t e ) 
and Ballistic (t < t e ). 

ballistic diffusive ergodic quantum 

t e t c t]{ =-^ t 


Figure 1. Schematic ordering of the diffusion time-scales. 

In metallic grains, in the ergodic window the spectral properties of the metals 
are universal, and described by the random matrix theory (universal conductance 
fluctuations), hence the second name of this regime. 

Let us try to make this description more quantitative. Let us introduce the 
retarded Green’s function G R (x, y\ E) =< x\(E — H + ie)~ l \y >. The time- 
Lourier transformation G R (x,y;t ) describes the amplitude, that the diffusing 
particle propagates from x to y in a time t, under the influence of the dynamics 
governed by some Hamiltonian H. The Hamiltonian is the microscopic source 
of the disorder. A classic example is the Anderson Hamiltonian [10] 

Ha = — tt - (V — ieA) 2 + V(x ) (27) 

Zm 

where < V(x) >= 0 and < V(x)V(y >= 6(x — y)/(2irp(Ep)t e ). Note 
that the disorder is static, i.e. time independent. We can now define a crucial 
concept of the return probability. This is simply the square of the amplitude 
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of returning to the same point x in a time t, averaged over the disorder. For a 
particle at a Fermi surface, it reads 

P(t) = lim-^ [ dXe iM (G R (x,y;E F + X/2)[G R (x,y, E F + A/2)]*) (28) 

The prefactor V comes from the translational invariance of the return probabil¬ 
ity, and the denominator 2itp(Ef) guarantees the normalization of the proba¬ 
bility to 1 . 

For a static random potential of the type above one can perform the averaging 
and then the integration, with the following result: 

P(x, y,t)=J2 Um ^°° _2__ e -| g-y\V±Dt (29) 

Hence 

p(t) = y e ~ DqH (30) 


or, equivalently, in Fourier space 


m = £ 

<? 


1 

—iX + Dq 2 


(31) 


In the above, while performing the averaging over the microscopic disorder, we 
integrated out fast degrees of freedom, getting the effective description in terms 
soft modes q,. Details of fast degrees of freedom arc now hidden in effective 
parameters like the diffusion constant D, and the obvious name "diffuson" for 
a soft modes q r is natural in the light of the form P(x. y. t ), being the Green’s 
function of a diffusion operator (dt + DV 2 ). 

Formally, for very large times (A —> 0), the return probability develops a pole 
(diffuson pole at q 2 = 0). However, for such times, quantum effects become 
relevant, introducing the natural cutoff 7 at the energy scale of the average 
quantum spacing A. This cutoff suppresses infinitely long diffusive orbits, and 
regulates the pole 

P(t) = e ~^ e~ DqH (32) 

q 


or, after Fourier transformation 

Pit) = Y — ---— T (33) 

“ — iX + 7 + Dq 2 

q 

In the next section, we will demonstrate, that all the above listed concepts arc 
directly applicable to QCD. 
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Euclidean QCD is diffusive in d=4 

Let us look at the Dirac operator in the background of some Euclidean gluonic 
configuration A 


+ im)q k = \ k {A)q k (34) 

as a "Hamiltonian" corresponding to the "eigenenergy" A. The imaginary con¬ 
stant shift im does not spoil the analogy. For static Hamiltonians in d = 1,2,3 
we know from the quantum mechanics, that the evolution in time t is governed 
by q k (t) = exp(iE k t)q k (0). Time t is dual to the energy E. We will parallel 
this construction here, introducing dual "time" r to the virtuality ("energy") A. 
Since the "Hamiltonian" Up(A) +im is by definition independent on this "time" 
{static in r), we can write down the eigenmode evolution in this "time" 

q k (r) = e^ k+i ^ T q k { 0) (35) 

corresponding to 4 + 1 dimensional Schroedinger-like equation with static (r 
independent) potential 

id T q(r) = ( UffA) +im)q{r). (36) 

Note that this construction does not modify any properties of the Dirac operator, 
and is basically equivalent to introducing the Schwinger proper time r. 

We will prove now, following [14], that the dynamics of the Euclidean d = 4 
Dirac operator in Schwinger time r parallels the dynamics of the usual diffusion 
d = 1,2, 3 in the real time t. This means, that we will identify, in the spectrum 
of the Dirac operator, four distinct regimes corresponding to quantum, ergodic, 
diffusive and ballistic regimes of the mesoscopic physics. 

From now on, we follow, step by step, the construction outlined in the pre¬ 
vious chapter. Instead of Green’s function for a diffusing electron we take a 
Green’s function for a quark in some unknown gluonic background, which is 
constant in Schwinger time r 

G R (x, y: E) —■> S F (x, y; A) (37) 

We can now define by analogy to (28) the return probability in time r for a 
quark at the zero virtuality surface 

P(t) = V , . / d\e iXr (TV S F St) (38) 

1 7 27 rp( 0 ) J \ f /qcd 

Here Sj? = Sp(x, y. m — iA/2), trace appeal - due to the y f , matrix structure of 
the Dirac operator and the averaging is done over the full QCD measure (21), 
representing the analog of Anderson Hamiltonian representing the disorder of 
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the system. The parameter m appeal's here due to the uniform shift im of the 
original spectrum of the massless Dirac operator. 

Till now the analogy was exact. But now we have to perform the averag¬ 
ing over the a priori unknown measure of the QCD. The first alternative is to 
choose some model of the disorder. Actually the analytical instanton model of 
Diakonov and Petrov [11] and the numerical simulation by Shuryak [12] are 
the first realizations of this scenario. We will return to this point later. 

Surprisingly, we can calculate P(r) without assuming any model of the 
disorder, but making use of almost exact low energy theorems of the QCD, 
introduced in Part 1. First, let us note, that due to the chiral properties of the 
Dirac operator 

< TV S F S\r >= - < TrS F (x, y; z)'y b S F (y, x; z* 75 > (39) 


where z = m — A/2. We may introduce also the isospin sources. One then 
immediately recognizes the similarity to the pion correlation function structure 


C£(x,y\m) = 






-E 

V Z - 


,iQ(x-y) 


Qt , 


<qq> 

F L 


Q 2 + Ml 


(40) 


where we used the pion dominance formulae (16). 

Second, we use another low-energy theorem (18), to replace the M n on 
the r.h.s. of the above equation by the current quark mass m. Third, we 
analytically continue the C u i b (x. y;m), by replacing the mass m by z = m — 
iX/2. We recognize then that the integrand in the return probability involves 
the analytically continued pion correlation function, so averaging is equivalent 
to the analytical continuation of the r.h.s. of C n 


P{x,y,r) 


d\e lXr C n (x,y; z(X)) 


2np 

]T f dXe iXT e iQ{x ~ y) 


F 2 


Q, 


l^Q L + 2m + iX 


(41) 


Finally, after integrating over A by the residue method, we see already the 
diffusive structure of P(x, y, r), and the return probability (limit x —> y) reads 

P{t) = e~ 2mT J2 e~ DQ2r (42) 

Qn 

or, after Fourier transforming 

f(A) = E 

Q[l 


iX + 2m + DQ 2 


(43) 
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We constructed in this way the precise analog of the diffusive return probability, 
with diffusion constant D = F 2 /\ < qq > | ~ 0.22 fm and the slow "diffuson 
modes" Q jt = Inn^/L, with p = 1,2, 3,4 and n /( integers. 

In fact, using the low energy theorems and the above construction, we im¬ 
plicitly integrated out the fast (gluonic) degrees of freedom, getting in this way 
the diffusion constant. The remaining dynamics of slow modes (diffusons) is 
the dynamics of the pions - longest wave excitations of the QCD. The Gold- 
stone nature of the pion manifest itself as a pole in the chiral limit and A —> 0 
limit. We see easily that the very heavy quarks do not diffuse at all (due to the 
exponential damping). Note also that for light but massive quarks we can read 
GOR as an expression of the coherence length of the QCD vacuum. Indeed, by 
definition, coherence length is 

L C oh = \J D/'y = 1 /M n (44) 

where we used that for QCD 7 = 2m and the expression for the diffusion 
constant calculated above. Coherence length of the QCD vacuum is related to 
the pion mass, so in the massless limit, pions are indeed much more vacuum 
modes than the pairs of bounded constituent quark and antiquark. 

Now we can identify the different regimes in the spectrum of the Dirac oper¬ 
ator. Since the Schwinger time r is dual to the eigenvalue of the Dirac operator, 
we immediately identify the analog of Thouless energy, here the Thouless vir¬ 
tuality [14, 15] 


A c = D/L 2 = F 2 / (T,Vv) (45) 

where | < qq > \ = S. Therefore the eigenvalues smaller than A c are ex¬ 
pected to belong to the ergodic (universal) regime, whereas larger than A c to 
the diffusive. Since the Banks Casher relation gives us an estimate of the mean 
level spacing, we know also the borderline of the quantum and ergodic regime. 
Finally, for very short times (large eigenvalues) the concept of the diffusion be¬ 
comes meaningless. It is not very difficult to argue, that the borderline between 
the diffusive and ballistic regime is given by twice the mass of the constituent 
quark, basically the mass of the p meson. Indeed, for times shorter than the 
required to travel one mean-free path between lumps of the gauge field the 
concept of dressing ( via multiple scatterings) of the current mass is void. 

The above hierarchy of scales could be summarized again by the cartoon, 
where we also remind about the thermodynamical ordering of the scales. On 
the basis of the analogy with the condensed matter systems, we may suspect, 
that similarly to the universal conductance fluctuation in the ergodic regime, 
we may see universal spectral fluctuations in the ergodic window of the QCD, 
hopefully described by some sort of random matrix theory. In the next chapters, 
we will see that this indeed is the case. 
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quantum ergodic diffusive ballistic 

A ~ y A c r}-t= 2 M q ~ >\ 

Figure 2. Disorder regimes in the eigenvalue spectrum of the Dirac operator. Note the 
thermodynamical ordering. 


Ergodic ( universal) regime of the QCD 

Let us consider in more detail the ergodic regime. We expand the sum 
over the diffuson (pion) modes valid in diffusive regime, i.e. we sum over the 
quadruples of integers Q M = (n\, ri 2 , n^, n^-jr/L 

P(t) = e“ 2mT (l + 8e- 4 ^ D/i2r + ...) (46) 

where the first term comes from all n t = 0 , second from all combinations of n 
of the type (1, 0,0,0) etc. Since D / L 2 = t c 1 . the return probability reads 

P ( t ) = e- 2mr (l + 8e- 47rV / Tc + ....) (47) 

For the times larger than t c all terms except of the first one vanish exponentially, 
so we obtain, that in the ergodic regime the return probability is simply 

P(t) = e~ 2mT (48) 

approaching constant for very large r. Only the softest modes (zero modes 
n\ = n ,2 = ns = 714 = 0 ) determine the return probability in the ergodic 
regime. This points that the properties of the ergodic regime are universal - they 
are independent on the details of space-time interactions, since the Goldstone 
bosons interaction involves derivative terms, which vanish in long range limit. 

This result was known since long ago, although rephrased in a different way. 
In the usual chiral perturbation theory, pion momenta arc of order 1/L, the mass 
of the pion scales therefore like 1/L 2 , and since m ~ Af 2 , the combination 
m 2 V ~ 0(1) is fixed. The systematic expansion based on this counting is 
the clue of the Weinberg [13] chiral perturbation expansion and was practically 
realized by Gasser and Leutwyler [16]. What is less known, that Gasser and 
Leutwyler looked also at finite volume, and what happens, when one keeps only 
zero modes of the pion propagator. Since the propagator in the final volume 
reads 
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the counting based on the concept of having mV fixed [17] (equivalent to 
M 2 V 0(1)) kills all the terms in the primed sum (non-zero modes) and 
leaves the zero mode. This is very different to the counting based on before- 
mentioned m 2 V = 1 principle (chiral perturbation theory). 

Gasser and Leutwyler managed to re-sum the contribution of all the zero 
modes based on mV ~ 0(1) counting, and obtained the exact partition function 

z(m) = f dUe VZTrm(U+Ul) ( 50) 

with U = exp(?'7rr) and S = \ < qq > |. This formula is universal: it depends 
only on the ways how chiral symmetry is to be broken spontaneously (choice 
of the measure dU of the Goldstone modes) and how chiral symmetry is being 
broken explicitly, by the + (Nf, Nj) representation in the exponent. 

In the ergodic regime any theory/model sharing the same global symmetries 
as QCD belongs to the same universality class, i.e. any theory/model leading 
to the return probability P{t) = e~ 2mr will exhibit the same universal spectral 
properties in the ergodic window. Since the information about the dimension¬ 
ality is being lost in the ergodic regime, this could be even a zero-dimensional 
version of the field theory. But field-theory in zero-dimensions is a matrix 
model. So in the ergodic window QCD is equivalent to a certain matrix model. 
Before we will unravel the details of this model, let us remind once more that 
this equivalence between QCD and Random Matrix Model (RMM) happened 
only for eigenvalues below A c . Since the edges of the ergodic window scale 
with volume as A ~ 1/V ( Banks-Casher) and as A c ~ 1/y/V ( Thouless 
virtuality in four dimensions) the ergodic, universal window shrinks to a point 
in the infinite volume limit. 

Where does the "color dirt" come from? 

At this moment we may start to worry, what plays the role of the diffusive 
dirt in the case of the QCD? QCD is a fundamental theory, so the "dirt" has to 
be an immanent feature of the QCD itself. Whereas we do not have an exact 
answer what is the color dirt in QCD, let us observe, that most of the localized. 
Euclidean (i.e. static from the point of view of Schwinger time) gluonic con¬ 
figurations do the job. Typical and perhaps the most natural are instantons [18]. 
First, they are the classical, localized stable Euclidean solutions of the QCD. 
Second, Dirac equations in the presence of instanton background possesses 
the chiral zero modes, therefore the instanton vacuum immediately provides 
a microscopic "hopping" mechanism from one instanton field to another and 
guarantees flipping the chirality at each "scattering" in "time" r. Quantita¬ 
tively, each instanton provides a seed of non-conservation of the chiral charge, 
by integrated form of the anomaly: 

5Q a = 2 N f Q 


(51) 
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where Q is the topological charge of the instanton. The instanton model in¬ 
volves basically two parameters, the average size of the instanton (p ~ 1/3 fm) 
and the concentration n ~ fm~ 4 (the typical density of instantons). Therefore 
every dimension-full quantity in the instanton model depends parametrically on 
the combination of the two scales. In particular', the smallness of the diffusion 
constant calculated in the previous chapter, finds a natural explanation as a di¬ 
luteness of the instanton medium. Note that the finite value of the condensation 
requires that thermodynamical limit V —> oo, N —> oo, with N/V fixed, where 
N is the number of instantons and antiinstantons. 

Each instanton vacuum configuration is a snapshot in a time r. Each snapshot 
violates the Lorentz invariance ("static", particular distribution of instantons and 
antiinstantons in four volume) and gauge invariance (each instanton freezes a 
direction in color). But the averaging over collective coordinates of the instan¬ 
tons (here over the centers of instantons and over their color measure) restores 
the gauge and Lorentz invariance. We would like to stress, that instantons are 
sufficient, but not necessary configurations to realize the diffusive scenario. 
Several other models may also provide chiral disorder, e.g. family of stochastic 
vacuum models [19]. Each of these models comes with a certain correlation 
length scale, which corresponds quantitatively to the diffusion constant. 

In other ways, the spontaneous breakdown of the chiral symmetry in QCD is 
a very robust phenomenon, comparing to the confinement. Usually, the models 
of color dirt either ignore confinement, or introduce it by hand, or as instantons, 
seem to be (at least naively) not related to confinement at all. Mysteriously, 
lattice evidence suggests strongly that both phenomena (confinement and chiral 
symmetry breakdown) are correlated and vanish at the same temperature. Since 
the confining configurations are the topic of vivid speculations, the fundamental 
understanding of the "color dirt" and chiral disorder is still missing. 

Part 3 

I analyze in more detail two regimes of the QCD, the ergodic one and the 
diffusive one. Then I confront few sample predictions with the "experimental 
data" obtained from lattice simulations. 

Random matrices - field theory in 0 dimension 

As stated before, in the ergodic regime any model obeying the global sym¬ 
metries of the QCD belongs to same universality class as QCD. Since this is the 
regime where only the constant pionic modes matter, we can ignore space-time 
dependence and stick to the field theory in zero dimensions, i.e. the theory 
where fields are numbers and do not have any space-time dependence. The 
exact form of such theory in the QCD can be inferred in numerous ways. Here 
we follow the historical route [11, 12, 20]. Imagine that we have two lumps of 
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the "dirt", e.g. instanton and anti-instanton, separated at the very large distance. 
Then, the Dirac equation in the field of the instanton has an exact fermionic 
right-handed zero mode cj)R, and the Dirac equation for the anti-instanton has 
a similar, but left-handed zero mode 4>l. When we decrease the distance be¬ 
tween instantons, the degenerated pair of zero modes is replaced by the pair 
of eigenvalues (T, —T), where the overlap T = f d 4 x(J R i<J)(l) /, depends on the 
distance and the mutual orientation of the instantons. Let us now add more and 
more instantons into this medium, and still work in a dilute gas approximation. 
The infinite fermion determinant, when calculated in the basis of left and right 
handed quark zero modes, is now approximated by the matrix of overlaps be¬ 
tween the I-th instanton and J-th anti-instanton T/j. The off-diagonal block 
structure comes from the chirality flipping mechanism. The diagonal blocks 
arc zero in the chiral limit. In this way, instanton picture trades the a priori 
unknown QCD measure (21) into the approximate measure 


<(...) >= / ..)e 

J N f 


0 

Tt(O0 


T(fii) 

0 


(52) 


where D* is the set of collective coordinates of the instantons ( positions and 
color orientations), and 5' g i ue is the gluonic paid of the QCD action saturated 
with the initial instanton vacuum ansatz. Let us truncate now all the space time 
dependence in the above action. We arc left with the model of the type 


Z = 


dTe - NaTrTT i det 


0 

Tt 


T 

0 


Nf 


(53) 


Here T is an A' by A' matrix built out of complex numbers (one could generalize 
the matrices to the rectangular ones as well). The gluonic measure is replaced 
by the polynomial measure with some potential v{T^ T), here being the simplest 
- the harmonic potential. The factor N in front of the potential guarantees, that 
each integration over T is appropriately weighted. The model has one scale, i.e. 
the "width" a of the Gaussian matrix measure. The partition function Z defines 
the chiral Gaussian random matrix model [21]. The name chiral comes from 
the off-diagonal block structure, which mimics the original chiral structure of 
the Dirac operator, \ljX A), 7 . 5 ]+ = 0. 

We will show now, that such defined model is exactly solvable for any finite 
N. The main problem is the integration over the measure (IT. We will use here 
the trick, similar to the change of the variables from Cartesian coordinates to 
the spherical ones. The simplest way to change usual coordinate variables is to 
look at the infinitesimal interval 

ds 2 = dx 2 + dy 2 + dz 2 = dr 2 + r 2 d0 2 + r 2 sin 2 6d(j) 2 (54) 

so the metric tensor is gif. = diag(l. r 2 , r 2 sin 2 0). Hence the Jacobian (the 
square root of the determinant of the metric tensor) reads J = \J det gif. = 
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y/g = r 2 sin 0. We observe now that the integrand of the partition function 
(53) depends only on the "radial" combination X = T^T, and we follow 
the trivial example above. Every her mi t!an matrix can be diagonalized by a 
unitary transformation. Since X is hermitian, introducing X = URW, where 
the unitary matrix U = exp (iH), H hermitian, and Rij = is positive 
diagonal, we calculate 

tr (dX 2 ) = tr(l7 t (dR + i[R, dH])U) 2 = tr (dR + i[R, dH]) 2 

= X! dr l + ~ r of dlli/r (55) 

k i^j 

Hence the metric tensor reads 


gik = diagQ^-^T, ..., (n - rj) 2 ,.... ) (56) 

tV N(N- 1) 

so det g tk = - rj) 2 and Jacobian J = = Yli<j(n - rj) 2 = 

A (R) 2 . Since the integrand and the Jacobian do not depend on the angles 
parameterizing the unitary matrices, we arc left with N integrations over the 
diagonal eigenvalues ("radial" variables) r, 

Z = j Y[driA{R) 2 ex.p(-Na^2ri)Y[r^ f (57) 


Let us note, that the expression A(i2) can be rewritten as a determinant (Van¬ 
dermonde determinant). 


A(i?) = n( r * “ r i) 

i<j 


1 r'i r\ 
1 V 2 r 2 

1 r N r% 


r 

r 


N -1 
1 

N -1 
2 


r 


N -1 
Af 


(58) 


[Exercise . Prove (58).] Since we can add to each row of the determinant 
an arbitrary combination of the other rows (without changing the value of the 
determinant), we can replace the original Vandermondian by the determinant 
build out of polynomials. 


A (R) = 

Poin) 

Po(r 2 ) 

Pi(r i) 
Pi(r 2 ) 

P-iin) . 
P 2 (r 2 ) .. 

.. Pat- 1 (ri) 

• • Pat-i (r 2 ) 


Po(r N ) 

Pi{r N ) 

P2 { r N ) • 

■ ■ Pw-i (tn) 


These polynomials arc a priori arbitrary, modulo the fact that the coefficient at 
the highest power is always 1. We choose now the polynomials in such a way. 
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that they arc orthogonal with respect to the measure. In our case we require, 
that 

/»oo 

/ dre- Nar r N fPi{r)P k (r)=5 ik (60) 

J O 

It is obvious, that in the case of Nj = 0 such polynomials arc Laguerre poly¬ 
nomials (modulo trivial rescaling of the coefficients, so the highest power is 
multiplied by 1) and in case of Nf / 0 arc associated Laguerre polynomials 

L k 1 . Due to the orthogonality of the polynomials all integrations in the parti¬ 
tion function factorize, and we can write down arbitrary (--correlation function 
between the k eigenvalues. Integrating over all eigenvalues except of the last 
one, we arrive at the average spectral density 


P(A) 


^ k\ 

t 0 n Nf + k+i) 


L { , Nf) L[ Nf) z N f + VV* 


N\z N f + h~ z r 
T(N f + N ) 


T (Nf) t (N f +1) _ AN f ) (N f +1) 


(61) 


where we used (in the second line) the Christoffel-Darboux summation formula 
and where z = NX 2 a is an argument of the polynomials, with A being the 
original eigenvalue of the matrix T. 

Let us now come back to physics. In the previous section we stressed, that 
the ergodic regime corresponds to the limit, where Vrn is fixed (in units of 
condensate £). This suggests, that the relevant limit for ergodic QCD follows 
from the chiral random matrix model with N A fixed (in units of the width a ) 1 . 
Technically this is achieved by so called unfolding procedure. We define the 
microscopic spectral density /Viicro(•“’') = lmi/v^oc jj-^f>(X = jP)- Using the 
asymptotic property of Laguerre polynomials 

lim -^L { %\x/N) = X - a / 2 J a {2^c) (62) 

iV—KX) iv 


we rewrite (61) in the microscopic limit in terms of the following combination 
of the Bessel functions 

Pmicro(s) = \v 2 s{J 2 Nf (crs) - J Nf -i(as)J Nf+1 (crs )) (63) 

This formula [23] represents the universal fluctuations ( universality class of 
chiral Gaussian unitary ensemble) of the spectral density in the QCD. The 
parameter a is fixed by tuning the dimensionless quantity N\a to dimen¬ 
sionless quantity kmS. It is important to mention, that this formula is uni¬ 
versal, i.e. similar for any polynomial potential [24] v(T'T) of the ran¬ 
dom matrix model. If we would not impose the limit NX = 0(1) while 
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N —> oo, instead of p m icro( s ) given by (63) we would have obtained the non- 
universal spectral density of the Gaussian random models, i.e. Wigner semi¬ 
circle p( A) ~ v/4 — A 2 . Colloquially speaking, all the universal fluctuations 



Figure 3. Microscopic spectral density for Nf = 0,1, 2 (solid, dashed dotted lines, respec¬ 
tively). 

are hidden at a single point of the Wigner semicircle, corresponding to point 
with coordinates (A = 0, p(0)). The "zooming" of this micro-dot by factor N 
unravels the fluctuating, universal pattern of microscopic spectral density of the 
chiral ensemble (Fig. 3). 


Crossover to the diffusive regime - chiral perturbation 


The universal (ergodic) regime ends at the scale corresponding to the Thou- 
less virtuality. For the scales larger than Thouless virtuality neglecting non-zero 
modes is no longer justified - all modes have to be taken into account. Let us 
look at the scaling properties of the eigenvalues. The first eigenvalue which 
falls off from the ergodic window (i.e. by definition, Thouless virtuality) has 
to scale [14, 15] (in units of A) as 


Ac D Fl FlV 2 2 

A A L 2 A| < qq > | L 2 ~ L 2 


(64) 


where we used the Banks-Casher relation and the definitions of the Thouless 
virtuality and the diffusion constant. 

It is remarkable, that the analogy with the diffusion can lead to a quantitative 
description of the spectral correlations in the diffusive regime of the QCD as 
well. Let us define two-point correlation density of states 


R(s) 


1 

< p( A) > 2 


<p(A - e/2)p(A + e/2)> - 1 


(65) 
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where s = e/A and, as usual, A ~ 1 /pV. Then we define the spectral 
formfactor 

K(t) = [ deR(e)e~ iet (66) 

27T ./ 

For times shorter than the Heisenberg time tu ~ A -1 , the standard semi- 
classical argument developed for mesoscopic systems by [25] allows to relate 
spectral formfactor to the return probability P(t ) 

K(t) = ^tA 2 P(t) (67) 

By definition, the two-level spectral correlation function, integrated from 0 to 
A measures the fluctuation X 2 (A) of the number of levels N (A) in a strip of a 
width A 


r<N> r<N> 

£ 2 (A) = / / dAidA 2 i?(Ai - A 2 ) (68) 

Jo Jo 

In the ergodic regime (for t —► oo) P(t) is a constant, hence K(t) ~ t. The 
variance X 2 grows as logarithm of A, as expected in the random matrix theory. 
However, for times smaller than t,i (eigenvalues greater then Thouless energy), 
the return probability is given by the classical diffusion result 


Pit) 


V 

(47T-Df) rf / 2 


(69) 


Then K(t ) ~ t 1 r/,/2 , which cotTesponds, after Fourier transforming, to 


R(s) ~ 


-Re A 2 ^ 
Q 


1 

(aA + iDQ 2 ) 2 


(70) 


For mesoscopic systems (e.g. tor disordered electrons in d = 3 metallic grains), 
this is the seminal result obtained by Altshuler and Shklovskii [26], demonstrat¬ 
ing, that the two-point correlation function comes diagrammatically from the 
two-diffuson exchange. The variance T?(E) reads then 



The second equation comes from definition of the Thouless scale and from 
E ~ f -1 = I) / L 2 (E). The power behavior reflects the fact, that in the 
diffusive regime diffusion of a particle with energy E is non-homogeneous, 
and takes place independently in the number of sub-blocks ( L/L(E)) d of the 
original block V = L d . 
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The above mentioned crossover from universal regime (logarithmic behavior 
of the variance) to diffusive regime (power behavior (E/E c ) 3 / 2 ) was confirmed 
by [27] in numerical simulation of the metallic regime for Anderson model with 
20 3 sites. 

On the basis of the diffusive scenario presented in these notes, we expect 
similar behavior for QCD, but with d = 4 and t —> r. In our case, the square 
denominator in (70) follows from the exchange of two pions in the double ring 
diagram corresponding to density-density correlation function, reflecting on 
disconnected quark susceptibility. Note that we have obtained this result on 
the basic of spectral analysis, and not by using a standard diagrammatics of 
the chiral perturbation. The diffusive regime of Altschuler-Shklovskii in meso¬ 
scopic systems is equivalent to chiral perturbation theory. The effective models 
of diffusons (sigma models) correspond to sigma models of the pions, and su¬ 
persymmetric formulation of the diffuson models by Efetov [28] resembles the 
family of so-called (partially) quenched chiral perturbation theories [29]. Finite 
volume QCD plays a role of a 4-dimensional quantum dot. 

For QCD, the diffusive regime should show the power behavior £ 2 (A) ~ 
(A/A c ) 2 , since the diffusion takes place in d = 4 ([14]). Recent lattice simula¬ 
tion [30] of the SU (3) QCD has confirmed this scenario, demonstrating clearly 
the expected crossover from ergodic regime (logarithmic behavior) to diffusive 
regime (power-like behavior) predicted for the QCD. It is illuminating to com¬ 
pare the character of the crossover in metallic grains (e.g. Fig. 3 in [27]) to 
crossover regime of the disconnected quark susceptibility (e.g. Fig. 8 in [30]). 

Fet me mention another intriguing analogy. In the case of Hamiltonians 
invariant under the time reflection, one does not differentiate between the di¬ 
rection, in which the loop contributing to return probability is traversed. Hence 
the naive contribution to the return probability is doubled, which corresponds to 
additional diffuson-like contribution in the formula for return probability. This 
quantum interference of the identical orbits traversed in opposite directions is 
called in condensed matter coherent backscattering (weak localization), and 
the collective diffuson like-excitation corresponding to this effect is called a 
cooperon, since the charges add during this interference. This phenomenon has 
also an analogy in QCD [14, 31]. In case of a real color group (e.g. SU (2)), 
we do not distinguish between the quark and antiquark, and on top of the usual 
pions (diffusons) we can form the collective states corresponding to qq pairs, 
so called baryonic pions - hence QCD-cooperons. The above mentioned lattice 
analysis [30] clearly identified the baryonic pions in quenched lattice simula¬ 
tion, in agreement with theoretical predictions based on the diffusive scenario 
of the QCD. 

In the light of the above results, it is tempting to speculate, that the ideas 
borrowed from mesoscopic systems may have much broader domain of appli¬ 
cations in the Euclidean QCD. Some preliminary analysis of the chiral disorder 
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influenced by external sources seems to show, that this is indeed the case [31]. 
The effects on disorder of quark chemical potential resemble a complex electric 
Aharonov-Bohm effect, breaking particle-antiparticle symmetry, accumulating 
two flux lines and leading to rupture of "baryonic" quark-antiquark pair at 
H = M,r/2. This is observed in quenched lattice simulations and in random 
matrix models [32], The phenomenon of persistence currents in disordered 
media finds an analogy to diffusion of light quarks in the presence of several 
Abelian Aharonov-Bohm fluxes. Some magnetic properties of chiral conden¬ 
sate can be explained by replacing the diffuson/cooperon trajectories by the 
four-dimensional Landau orbits. Low energy theorems obtained in [33] can 
be interpreted as negative magnetoresistance of the QCD vacuum (quark con¬ 
densate grows with the magnetic field). Low temperatures correspond to the 
replacement of initial cube L 4 by asymmetric box /)L'\ and the diffuson modes 
reproduce the lowest temperature corrections for the chiral observables. 

An extremely interesting case corresponds to the situation, when the critical 
temperature is reached. In the case of a second order phase transition, pion 
wavefunction, susceptibility and the condensate undergo the following scaling 


= £^ //3 , ** = S 1 " 5 , £ = m 1/s (72) 


where I used the standard notation for the critical exponents. These effects 
have an obvious quantitative effects on the diffuson (pion), modifying the re¬ 
turn probability. In particular, the probability of return at the critical point 
tends to universal behavior P(t) = (th/t) 1 ^ . This has to be contrasted with 
the vacuum result P(r) —» const, in the ergodic regime and vacuum result 
P(t) —> (r c /r) 2 in the diffusive regime. Finally, let me point that the con¬ 
cept of the return probability at the Schwinger time allows to check scenarios 
beyond the rigorous theory of phase transitions. For metal-insulator type tran¬ 
sition, P(t) ~ r 1 The multifractal exponent // estimated on the basis 
of chiral disorder is rj /4 = 0.057 [31]. The smallness of 7] makes the return 
probability to look "diffusive" with d = 2. Finally, the asymmetry of the box 
may cause the appearance of two diffusion processes: "temporal" one and the 
"spatial" one. Since the spatial diffusion constant D s is small in the high tem¬ 
perature phase (typically, D s /D t ~ 0.1 at T ~ 180 MeV), the asymmetry in 
conduction properties may cause a "percolation" from d = 4 to d = 1, with 
one-dimensional diffusive behavior P(t) ~ 1 /\Jr. These speculations are 
interesting from the point of the view of lattice spectral analysis in the vicinity 
of the critical temperature, since the return probability in the Schwinger time 
can be explicitly expressed in terms of the eigenfunctions and eigenvalues of 
the Dirac operator. 
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Threefold ways and lattice spectra 

Let me finish these lectures reminding the classification of the global sym¬ 
metries of the Dirac operator, hence the classification of the chiral universality 
classes. 

For complex quarks, the Dirac operator has no symmetries. For real quarks 
(two colors), Dirac has an anti-unitary symmetry, [ 727402 K, Upp] = 0, where 
K denotes complex conjugation. For adjoint quarks Dirac operator has another 
anti-unitary symmetry, [7274 K, HJ)a\ = 0, where subscripts F. A denote the 
fundamental and adjoint representations, respectively. This completes the list, 
forming the threefold way, how global symmetries can be realized for the Dirac 
Hamiltonian [34]: 

I. For complex quarks with Nf flavors, the pattern of spontaneous breakdown 
of the chiral symmetry corresponds to SU(Nf) x SU(Nf) —► ,3’(/(Ay), and 
the random matrix realization is obtained by filling the matrix with complex 
numbers (chiral (denoted by x) Gaussian unitary ensemble (yC'ue)). This is 
the case of the QCD, analyzed in these lectures. 

II. For real quarks with Ay flavors, the pattern of spontaneous breakdown 
of the chiral symmetry corresponds to SU(2Nf ) —> Sp(2Nf), and the random 
matrix realization is obtained by filling the matrix with real numbers ( chiral 
Gaussian orthogonal ensemble ( xGoe )). This is the case of QCD with two 
colors. 

III. For adjoint quarks with even Ay, the pattern of spontaneous breakdown 
of the chiral symmetry corresponds to .S'(/(Ay) —r 0(2Nj), and the random 
matrix realization is obtained by filling the matrix with quaternions ( chiral 
Gaussian symplectic ensemble ( xGspe )). This is the case of QCD with adjoint 
fermions, for any number of colors. 

For each of these ensembles, one can obtain exact spectral formulae. Histor¬ 
ically, the first one was (63), obtained in [23]. Technically, the complex case is 
the easiest one. 

This triad is distorted on the lattice, due to the known problems of incom¬ 
patibility of having chiral fermions on the discrete lattice with local action 2 . 

For example, the additional symmetries of Kogut-Susskind Hamiltonian 
cause, that for two colors the lattice universality class is xGspe, and not, as ex¬ 
pected from II, xGoe. Historically, the microscopic spectrum for this ensemble 
was the first one explicitly confirmed by lattice, in [35]. The formula (63) was 
confirmed successfully on the lattice (Nf = 0 case) only recently [36]. 

Nowadays, there exist an impressive list of theoretical predictions for several 
spectral (also higher-point) correlations and plethora of lattice evidence for 
various combinations of gauge groups, fermion representations, topological 
sectors, with quenched and unquenched and "double scaled" determinants etc. 
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We refer the curious reader to original literature. In all cases, the agreement 
with RMM is excellent. 

There is also a growing interest in moving toward the diffusive regime. First, 
the scaling (64) was confirmed by lattice study [37]. In the previous subsection 
we mentioned the measurements of the crossover between the ergodic and 
diffusive regime. They arc also first measurements on the chiral parameters, 
motivated by the spectral properties of the Dirac operator. 

Further reading 

An excellent review on chiral dynamics is [38], where the reader will find 
deeper justification of several statements made in Paid 1. The original arguments 
on diffusion in QCD , based on rather compact paper [14] and the sequels [31], 
will become more obvious after reading condensed matter literature on this 
subject [7, 8]. The consequences of Banks-Casher relation are discussed in [6]. 
Basic informations on the phenomenology of the instanton vacuum could be 
found in [18]. Mehta’ s book [39] is a classics in random matrix theory, for 
the reviews on the chiral ensembles we refer to [40, 41]. Readers interested in 
broader aspects of RMM we refer to [42, 43]. I did not discuss here at all the 
schematic use of random models in QCD phenomenology, concentrating here 
on exact, quantitative predictions in the ergodic window. Sample applications 
could be found in a review [44] and references therein. Finally, a guidebook 
to advanced details (supersymmetric technique, replicas etc. ) of the partially 
quenched perturbation techniques could be found in recent reviews [41, 45] 
and references therein. 

Conclusions 

In these lectures, I tried to demonstrate in a physical way how much the 
concepts of the spontaneous breakdown of the chiral symmetry share with the 
concepts of the disorder in mesoscopic systems. This should comply with the 
introductory character of these lectures. I perhaps omitted several important 
works to quote, and I apologize for this. The main idea of these lectures 
was to introduce a non-expert to the fascinating field of chiral disorder in the 
QCD. I therefore sacrificed several advanced technical details and mathematical 
aspects, which are often necessary to prove the analogies presented here in a 
strict way. Luckily, there exist reviews, some quoted above, which exhaustively 
can guide the reader through the technicalities of random matrices and various 
versions of chiral perturbation theories adapted for finite volume systems. They 
also provide the more complete bibliography. 

I pointed, that the ergodic (universal) domain, is to large extent understood 
and documented in an impressive way. Universality of the ergodic regime 
helped lattice QCD to check to what extend chiral properties (including the 



LECTURES ON CHIRAL DISORDER IN QCD 


27 


global features e.g. zero modes) of the quarks arc reflected by simulation. They 
also allowed to quantify the discrepancies, in the light of the exact analytical 
predictions of the RMM. I tried to emphasize, how challenging is to look at 
the spectra of the diffusive regime. The lattice analogs of the diffuson theories 
open a way to systematic extraction of the effective parameters of the chiral 
Lagrangian, following step by step the Weinberg expansion. (The first example 
is a pion decay constant). It is also possible to study the dependence of these 
constants of the chiral Lagrangian on the external parameters and it is important 
to see what happens to them when the critical regime is approached. It is 
also plausible, that the spectral analysis in the diffusive regime may serve as 
an important tool for studying the nonperturbative nature of the gluon fields 
and their quantitative rearrangements in the vicinity of the phase transition, 
hopefully shading some more light on the nature of the lumps of the "color 
dirt". 

A more general message, which the reader should infer from these lectures 
and from several other lectures at this school is, how broadly and how fruitfully 
the theory of strong interactions can borrow from the concepts and ideas of the 
condensed matter theory. 
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Notes 

1. The rigorous proof of the equivalence of RMM limit and chiral partition function is non-trivial [22]. 

2. The original classification can be however reproduced using nonlocal actions. 
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